Coherence and incoherence in extended broad band triplet interaction* 
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In the present analysis we study the transition from coherent to incoherent dynamics in a non- 
hnear triplet of broad band combs of waves. Expanding the analysis of previous works, this paper 
investigates what happens when the band of available modes is much larger than that of the initial 
narrower combs within which the nonlinear interaction is not subjected to selection rules involving 
wave momenta. Here selection rules are present and active, and we examine how and when coherence 
can be defined. 
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I. INTRODUCTION 

Wave triplet interactions model a vast number of cases 
where nonlinear wave dynamics of physical systems can 
be described in terms of three dominant modes. The in- 
teraction is seen in a variety of situations, ranging from 
three wave interactions in laser-plasma and optical sys- 
tems to pulsar emission of electromagnetic radiation, in- 
cluding wave interaction in fluids and in several other 
settings 

The conservative interaction, which will be our focus 
here, is more easily handled when the interaction involves 
only the three pure modes of the triplet. However, a more 
realistic view should allow for a microscopic description, 
where each of the pure modes is replaced with a comb 
with many submodes. This has been done in a number 
of papers [1, H, 0, Q where several results have been de- 
rived along recent years. The main lesson one learns is 
that the dynamics can be correctly described in terms of 
three single or central modes, as long as nonlinearities 
are strong enough to lock all submodes into a single co- 
herent mode. For practical purposes the coherent modes 
can then be viewed as the pure modes of the triplet inter- 
action. On the other hand, if locking is not effective, each 
of the submodes follows its own linear dynamics and co- 
herence is lost. Random phase approximations can then 
be invoked to analyze the problem Q, but the concept 
of a pure triplet has to be abandoned. 

A recent paper Q shows how combs of modes can be 
very naturally formed in a wave system: the essential 
requirement, as we shall review, is that the nonlinear in- 
teraction takes place under spatially inhomogeneous con- 
ditions. When the inhomogeneity is present, wave vector 
matching among the interacting modes needs not to be 
exact since it includes the reciprocal vectors of the in- 
homogeneities. What happens then is that even if the 
initial conditions involve a small number of wave modes. 
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in a very short time interval the initial modes scatter off 
the inhomogeneities, creating groups of many modes, the 
initial combs. Another but equivalent way to see how 
combs are related to inhomogeneities is to realize that in 
the interaction of tightly packed group of modes, neigh- 
boring wave vectors cannot be properly resolved in finite 
size spatial scales, a typical occurrence in experimental 
settings [l^, UH . In this case whole groups of modes with 
similar wave vectors are altogether excited forming the 
combs. The interaction acquires the aspect of a mean 
field theory, where modes of one comb interact with av- 
erages taken over modes of the remaining combs [T^j. 

In the past, models for wave combs were based on 
combs with fixed number of modes. Once the combs were 
formed, submodes could evolve in time, but always pre- 
serving a prefixed total number within each of the combs. 
A recent analysis shows that combs with prefixed 
number of modes cannot actually maintain this number 
if the interaction takes place in a homogeneous environ- 
ment. As one may conclude from the comments above, 
this happens in virtue of the fact that homogeneity is 
unable to create a natural wave vector scale which could 
accommodate a given finite number of modes. Ref. [l^ 
indeed shows that as the wave dynamics develops, more 
and more modes are gradually excited and included in 
the interaction. 

This leads us to the central question of the present 
analysis, namely can the wave interaction in inhomoge- 
neous settings be well described with combs of finite num- 
ber of modes? We shall see that the answer depends on 
the time scales and the wave vector scales one is inter- 
ested in. 

The plan of the paper is the following. In §2 we first 
define a convenient interaction model allowing for an in- 
homogeneous environment and explore how the model 
can be used to create the picture of interacting combs 
with fixed number of modes, simultaneously analyzing 
its inherent limitations. In §3 we examine what hap- 
pens when the constraint of a constant prefixed number 
of modes is relaxed. In S4 we summarize our results. 
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II. THE MODEL 

The investigation starts as we consider the set of fully 
dimensionless space time equations governing the decay 
of mode "i" into modes and "5": 

idtai{x, t) + ivgidxttiix, t) = s{x)a2{x, t)a3{x, t), (1) 
idta2ix,t) + ivg2dxa2{x,t) s{x)ai{x,t)a3{x,t)* , (2) 
idta3{x,t) + ivg3dxa3{x,t) ^ s{x)ai{x,t)a2{x,t)* . (3) 

Set (HI) - (O actually describes the slow modulational 
dynamics for the complex wave amplitudes ap{x,t) (p — 
1,2,3) of corresponding carrier waves whose frequencies 
and wave vectors are matched. The combs are thus 
the multitude of sideband modes forming around each 
of three high frequency carriers, = —1, and the real 
function s{x) is the spatially dependent form factor in- 
troducing inhomogeneity in the problem. Function s{x) 
could be typically associated with inhomogcneous density 
distributions in plasma systems for instance. 

Let us first of all see how the classical picture of combs 
with given number of modes can arise from the basic 



set. We first need a structure for the function s{x). We 
define it as an even function centered at x = and a 
characteristic half width Is, as follows: 

s = s{x/ls) = s{\x\/h); s{\x\/ls » 1) ^ 0, s(0) = 1, 

(4) 

where for mathematical convenience, and with no loss of 
generality, we assumed a scaling that renders s(0) = 1. 
This kind of function restricts the effective interaction 
region as commented in the Introduction and can be used 
to introduce the basic wave vector associated with the 
inhomogeneities of the system in the form kg ~ l/^s. 
Now we write each of the waves ap{x, t) as combs of many 
modes 



ap{x,t) = J ap{Kp,t)e''^'''^dKp, (5) 

where Kp denotes the wave vectors of submodes within 
each comb. 

Spatial Fourier analysis of set H]) - ([3]) produces the 
following group of equations for the various submodes: 



iai(Ki) = Wgi Kiai(Ki) + / s(ki - K2 - K3) a2('«2)a3(K3)c^K2 d'«3, (6) 

iap{K2) ^ Vg2 K2a2(K2) + / s (ki - K2 - K3) Ci (ki )a3 (K3) ^^3 , (7) 



iaaiKs) = Vg3 Ksa^ins) + s(ki - K2 - K3) ai('«i)a2(K2)*c^Ki d^a, (8) 

J Ki ,K2 



r 

with 



s{x/ls 



1 

2^ 



s{x/ls)e '"'^dx, 



(9) 



s{ks) also even. One thus sees from the second of Eqs. ([9]) 
that in general, wave vector mismatches of magnitudes 
up to \ki — K2 — K^lmax ~ T^/h amoug the interacting sub- 
modes are allowed. If one defines a band width A in the 
form —A/2 < Kp < A/2, one concludes that all modes 
initially placed within the bands will interact simultane- 
ously, with no constraints due to selection rules, provided 
tt/Is ~ 3A/2. We shall refer to this regime as the regime 
of democratic interaction because selection rules are not 
operative here; under this regime, any three modes within 
the bands are coupled with the same strength. If Zs ^ 00 



one recovers the matched selection rule ki = K2 + K3, 
but for finite l^'s any triple of modes within the bands 
are connected. The approximate dynamics of bands can 
be obtained if one assumes s{ki — K2 — K3) ~ s(0) for 
\Kp\ < A/2, discarding all modes outside the combs; we 
note that under this approximation, and considering the 
normalization choice s(0) = 1, the first of Eqs. (O in- 
forms us that s(0) ^ 1/(3 A). In this case, and mov- 
ing into the discrete version of our continuum equations 
with Kp=i,2.3 — > = 2TTm/L = rriKL ("to" is an integer 
denoting the modal number), dKpCip^Kp) — KLdp{Kp) — 
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{2tt / L)ap{Kp) cipm, and L as the system length, one 
arrives at the set aheady explored by various authors 



iaiq = VglKqaiq + — — } a2ma3n, (10) 
m.n 

ia2m = Vg2Hma2m + y^^^-lq^-ln^ (H) 

q,n 

q.in 

To obtain set HO]) - ([12]) the prefactor kl/(3A) = 
(1/3) (k^/A) of the nonlinear terms in the discrete ver- 
sion is written as (1/3) /(1/-/Va), N^^ = A/kl being 
therefore a measure of the number of modes composing 
the combs in the Fourier reciprocal space; the factor of 
3 can be absorbed into convenient rescalings. As men- 
tioned, set Uni) - IT^ comprises the classical form of 



the broad band triplet interaction, where selection rules 
among the wave vectors are absent in virtue of finite size 
of the interaction region. Several interesting results have 
been obtained, the most prominent of which concerning 
the competition between the linear and nonlinear terms. 
If the linear band width terms associated with the group 
velocities are absent, one shows that in steady state the 
wave systems oscillates with a single nonlinear frequency 
Q.. If |ri| is larger than the largest linear frequency WgA/2 
(when unnecessary, modal and comb subindexes are oc- 
casionally suppressed to simplify notation), a phase lock- 
ing mechanism is present, preventing an initially coherent 
comb to decohere. In general, when a linear band width 
is present a time propagator g{t) can be constructed for 
the total amplitude, or macroscopic field of each comb 

Ap = ^ dpj, (13) 

3 

in the form Q 



g{t) — l^M2 — «e™"*'c'K — —i — igt\ ^ (time domain) 

(14) 

giij) = '"[(--".A/2)2-ln[(^+..,A/2)^] _ ..Sign(.,A/2-^^)+.^.Sign(.,A/2+^) (f^g^^g^^y domain). 



The factor —i exp{ivg Kt) in the time domain expres- 
sion is essentially the propagator for the microscopic 
mode with wave vector k, and the total propagator is 
obtained through an integration over the whole comb. 
If in the second of Eqs. one identifies the Fourier 

frequency tu with the dominant nonlinear frequency 17, 
the conclusion is that a dissipative-like term arises when- 
ever \Q\ < VgA/2. In extreme nonlinear cases with 
> VgA/2, coherence is preserved. In fact, a rela- 
tively straightforward procedure involving expansion of 
g around and a Fourier inversion from frequency to 
time domain, allows to write a coupled set for the the 
macroscopic fields @ which gives a good qualitative view 
of the dynamics in the democratic regime: 

iAi W + ^^2^3 (15) 

tA2 « /32^2 + ^AiA* (16) 

lAs ^ + ^AiA*2, (17) 

where jS (wgA)^/(12 il) iivgA < fl, and (3 -iVgA if 
VgA ^ ri. One sees that given the autonomous aspect of 
set P3|) - (IT71) one predicts decay (shrinking of volumes 
in the corresponding phase space) if A becomes larger 
than the nonlinear frequency. 

We shall obtain explicitly for some cases, but let us 
first dwell on the role of the width A. It is a fixed quantity 



which corresponds to one third of the total interaction 
range defined by the form factor s{ks). The width A con- 
tains a number = A/k^ of modes which in the past 
were supposed to be the only active modes of the wave 
system. However, the traditional model set (fTO)) - (fT2|) is 
only an approximation to the full nonlinear system ^ - 
^ , where one deliberately discards all modes outside the 
comb of the given width A. The assumption looks right 
because, as mentioned, modes within the comb are ex- 
pected to be more strongly and more quickly excited than 
modes outside. However, when one looks at the full set of 
equations there is always a nonlinear coupling which may 
eventually interlace and excite all modes, even those not 
initially placed inside the combs. In a real system with 
a band extension naturally much larger than the width 
A, the propagator for the entire macroscopic field should 
be rewritten as in Eq. but with A replaced with 

At, the latter quantity representing the total band width 
available to the modes. Thus, even if \n\ > Vg A/2, one 
might still have jfij < Vg At/ 2, a situation where coher- 
ence decay might be present. Of course, if one takes At 
as the full band width, and At > A, not all modes will 
interact democratically and selection rules shall reappear. 
In that case, previous results must be re-evaluated. In 
particular, from the stand point of macroscopic modes, 
the systems ceases to be autonomous since the nonlin- 
ear terms can no longer be written only in terms of Ai, 
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A2 , and A3 . Therefore one cannot prove or disprove that 
volumes in the phase space of the macroscopic modes are 
shrinking, as it happens with the approximate form given 
by Eqs. (fTS]) - (fTT]) . Nevertheless a dissipative term is 
present and the macroscopic modes are likely to decay in 
time - this is what really happens as we show next. 

III. A MORE ACCURATE VIEW: THE 
EXTENDED BROAD BAND INTERACTION 

As said, the full set © - I® is equivalent to the its 
counterpart spatial set H]) - ([3]). The connection is rel- 
evant because if one discards space derivatives exact so- 
lutions can be obtained. These exact solutions form the 
basis for further progress as one includes the space deriva- 
tives. 



A. Neglecting space derivatives 

Taking Vgpdx in the Eqs. ([T]) - fS]) a station- 
ary solution can be obtained in the form Ap{x,t) — 
Pp{x)exjp{i4)p{x,t)), where 01 = -2pi{x) s{x)t, 02,3 = 
—pi{x)s{x)t, p2,3 — V2pi{x), and where pi{x) is an 
arbitrary x-dependent function; we note that in the sta- 
tionary state phases depend on time, but amplitudes do 
not. Once pi{x) is defined, the complete solution is au- 
tomatically found. And once the space time solution is 
found, Fourier transforms can be used to move into the 
reciprocal space. To further specify the system with ini- 
tial conditions localized both in the real and reciprocal 
spaces, we shall make the following choice for the combs 
and the form factor s in the spatial representation: 

pi{x) = po exp (-x^/^p) , 
six/Is) ^exp{-~xyil) . (18) 

The comb thus defined has width ^ 1//^ in the reciprocal 
space and in order that its modes interact democratically 



at least initially, we require 1/1^ > 5/lp as explained ear- 
lier - in all numerical work we actually take l/l^ — 3/lp. 
We shall also assume that At ^ A, and write for the 
exponential distribution A — 2iT/lp. 

Independently of the choices we make for pi{x) and 
s{x) we are already in position to define coherence in the 
present case. We simply note that since 

ap{x,t) =^apje'^''', (19) 
3 

it is true that the macroscopic fields Ap introduced earlier 
in Eq. (|13p obey simple expressions - we write down the 
one obeyed by Ai: 

^1 = ^aij = Poexp(-2ipoi) • (20) 
3 

From the equation above we see that the macroscopic 
field oscillates harmonically with frequency il = —2pQ 
and with constant amplitude po- This is what we shall 
take as a coherent state: a non-decaying macroscopic 
mode oscillating with constant amplitude and constant 
frequency. The question that poses itself here is to deter- 
mine how many microscopic modes actually participate 
in the coherent state. In other words, would be true to 
assume that only the modes inside the initially defined 
combs are active? At a first glance one might suspect the 
answer would be positive since those are the modes in- 
teracting more strongly in the system. However, we had 
already pointed out that due to the nonlinear cascading 
structure of the interaction, some energy may flow from 
low to high wave vectors; and this is what actually hap- 
pens. This can be seen more formally with help of some 
tools. Keeping focus on the first comb, one first chooses 
a range D defined by —D/2 < k < D/2 and performs a 
partial summation over its internal modes. 



D/2 



ai{Ki)dKi 



-D/2 



271 



D/2 



-D/2 J - 



^ Pi{x) exp [— 2ipi(x) s(a;) i] dxdni. 



(21) 



r 



For a finite band D, the integral over ki, performed 
firstly, yields a delta-like structure as a function of x, 
with height D/27r and width 2tt/D. If one supposes 1/D 
small, the remaining integration over x can be done with 
help of a saddle approximation near x = Q where the 
space derivative of fields and form factor vanish. The 



final result can be written in the form 
Erfi 



X, 



D 



-2i pq t 



P^D- 



g3W47r^2pots"(0)/2/i:> 



(22) 

where s" = d'^s/dx'^, where Erfi(x) denotes the imag- 
inary error function as a function of argument x, and 
where we recall that s{x) varies faster than pi{x). We 
see that all depends on the behavior of the imaginary er- 



ror function for large and small arguments. If |x| <C 1, 
Erfi(x) ~ X and if |x| ^ 1, Erfi(x) ~ i- One therefore 
concludes that 



Constant when t < l-^ D"^ / 2 tt^ po 



when t > 15^/2 ttVo- 



(23) 



In other words, given a range D there is an intrinsic lim- 
iting time for coherence, 



2 ;2 



TD = 



(24) 



where by intrinsic we understand the limiting time ob- 
tained in the absence of the linear frequency band width, 
i.e. by taking Vg = 0. We know from our discussion re- 
garding Eq. ^ that Is ~ tt/ (3 A/2), so, the intrinsic co- 
herence time for modes within the original packet would 
be given by T£)=a ~ 1/po which is relatively small since 
this is essentially the period of the nonlinear wave. Our 
conclusion is that the initial packet can be hardly called 
a coherent structure even in the absence of the frequency 
band width. The collection of modes that could be seen 
as a coherent structure is anyone where D ^ A. In that 
case it is still true that decay will be present, but for all 
practical purposes tu would be so much larger than the 
period of the nonlinear wave that a physical setting or 
equipment resolving modes up to k ~ I? would perceive 
the wave system as coherent. 

A second important time scale has to be defined for the 
wave system. It is the time scale of excitation of individ- 
ual modes in the reciprocal space. Looking again at the 
first comb - reasonings are similar for the other two, we 
first recall the expression Ai(x,t) = pi{x) e'>q) (i(l)i{x,t)) 
for the steady state field. pi{x) is constant in time and 
the phase 4>i{x) = —2pi{x)s{x)t depends both on the 
spatial coordinate and time. If one evaluates the phase 
gradient d(j)i jdx and look at the maximum of this quan- 
tity as the largest wave vector involved in the dynamics, 
one derives the relation 



(2poA.)i, 



(25) 



which shows that the packet spreads over the recipro- 
cal space at a rate ±2po/^s- The time for excitation of 
any particular wave vector Kmax is thus Te^j^cij^max) — 
K.maxh/'^PQ- If we take Kmax — D/2, wc sce that for 
the typical case > 1, td » Texci^max = D/2), 
which means that in the absence of linear frequency band 
widths, coherence time of a packet of range D is in gen- 
eral much longer than the time required to activate the 
modes at the borders of the packet. 

In Fig. [1] we display the contrasting behaviors for lo 
in the cases Dig ^ 1 and DZs ^ 1- In the simulations 
we integrate set (P) - ([3]) with a pseudo-spectral method, 
using a grid of length L = N ~ 2^^, N denoting the 
number of nodes which for scaling simplicity is equal to 
the length. In all numerical analysis we use po = 1] the 
choice is not restrictive because field scales can always 
be absorbed in space and time. Considering Ig = 2^ and 
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FIG. 1: Id as a function of time for D Is = 2 n x 2^ in (a) and 
for D Is — 2ttx2^ in(b). In (c) we show the times series for the 
real part of the borderline mode with wave vector k — D/2, 
Is = 2^ Ip = 3ls. The group velocity is zero for all waves and 
all quantities are dimensionless. 



Ip — 3 Is, panel (a) displays the case D = 2tt/L x 2^, 
for which Dls = 2tt x 2^, and tjj ^ 32, while in panel 
(h) Dls = 2n X 2^ for which td ^ 8200. Panel (a) re- 
veals a fast decay, but coherence is far more persistent 
in panel (b). It is noticeable that in panel (b) function 
Id, although initially laminar, develops slight modula- 
tions after a very sharp instant along the time axis. This 
very sharp instant corresponds to t = TexciD/2). Indeed, 
the excitation time reads Texc{D/2) = 100.5 in this par- 
ticular instance. This is confirmed in panel (c) where, 
for the same parameters of panel (b), we show the time 
evolution of the real part of mode with wave vector D/2. 
The superscript "r" means "real part" and the submodal 
index d reads d = {D /2)/{2tt/ L) = 2^^ in this case, as 
defined in the context of the discrete equations, Eqs. (fTU]) 
- (|12p . We emphasize that as mode k = D/2 is excited, 
coherence of the packet D, although undergoing a mod- 
ulational process, does not decay. 

Of course, the presence of a band width for the linear 
frequencies may change the entire picture, and this is the 
subject of the next section. 



B. The effects of space derivatives and the 
associated linear frequency band width 

Since the full nonlinear system is not autonomous from 
the perspective of macroscopic modes, one cannot make 
very formal predictions about coherence decay due to the 
frequency band widths, like we did in the approximations 
leading to Eqs. - P?)) . However, some estimates can 
still be made. 

Let us consider our expression (|25p for the maximal 
wave vector involved into the dynamics. When Umax 



reaches the value Hr corresponding to the resonant fre- 
quency, VgiKr = 2poj coherence is expected to be lost, 
but now due to resonant effects. Under this circum- 
stance, the largest excited linear frequency would become 
comparable to the nonlinear triplet frequency fl, and co- 
herent nonlinear effects would be no longer dominant. 
The time to attain resonance, let us call it = Texcii'^r), 
can be obtained as one uses K^ax = K,r in Eq. 

Tr = ls/Vgl. (26) 

A given collection of modes of range D will remain co- 
herent as long as t < Tcoh = niin{T£), r^}. To illustrate 
this point, let us take the case analyzed in the panel (b) 
of Fig. [TJ In that case to is large and we do not ex- 
pect to see coherence decay soon if the linear frequencies 
are absent. But now let us add a frequency band width 
with Vgi chosen such that a given mode of the spectrum 
becomes resonant with < ; we achieve this require- 
ment with Vgi = 1/(k/2^) which yields ~ 50 < t^. 
For completeness we take Wgi,2 = which corresponds 
to one wave moving relatively to the other two. The 
setting would be of relevance to Brillouin scattering, for 
instance, where two electromagnetic waves with the same 
group velocity interact with a slower ion wave; we would 
be examining the process in the frame where the electro- 
magnetic wave is stationary. The resulting dynamics is 
then displayed in Fig. [21 where one clearly sees a fast 
decay whereas for Vg = Q one sees persistent coherence 
as previously shown in Fig. [^b). 

Expressions (|24p and (|26p therefore provide us with a 
simple tool to make estimates on the circumstances al- 
lowing coherence to be seen in the nonlinear triplet sys- 
tem. Once one has defined an extended comb distributed 
over a range —D/2 < k < D/2 of wave vectors with 
D/ A ^ 1, and once one knows the group velocity Vg for 
this particular class of wave, the coherence time can be 
obtained. 



IV. FINAL REMARKS 

In this paper we developed a technique to investigate 
coherence in nonlinear triplets, when the available band 
of modes is much larger than that of the initial combs. If 
modes remain restricted to their initial combs, the series 
of approximations outlined in §2 allows to describe the 
system as an interaction of macroscopic modes. In the 
presently studied case, one cannot resort to these approx- 
imations because initially low amplitude, idle modes out- 
side the initial range will be gradually excited at a rate 
~ 2po/ls, whenever the whole available band is larger 
than A. Coherence in this, perhaps, more realistic case 
is a little more involved subject to define. One first de- 
fines the range D of interest. The range has an intrinsic 
coherence time td defined in the absence of any frequency 



b; 
tl 

ta 




time 

FIG. 2: Coherence decay due to the resonant effect. Param- 
eters are those of panel (b) of Fig. [1] with exception of Vgi 
which here reads Vgi = 1/(2^ kl), defining a resonant time 

Tr « 50. 



the range D is defined, one has to look at the excitation 
time Tr of the resonant mode, which does not necessarily 
belongs to the range D; we found that Tr = h/vg. Gath- 
ering together both time scales, the final conclusion is 
that the coherence time Tcoh satisfies Tcoh = min {td, Tr}. 
We have also observed and stressed that coherence gains 
some substantial meaning only when several nonlinear 
oscillations occur prior to to- Since in our normalized 
variables the period of the nonlinear oscillation is ~ f /poj 
one concludes that the dynamics resembles a nonlinear 
phase locking process only when £) 3> A and Vg A <^ p^. 

Let us connect our results with those of previous works. 
Our macroscopic model does not look into fine micro- 
scopic scales of size, say Imic, where discrete effects be- 
come relevant. Therefore an upper limit Dmax ^/Imic 
does exist beyond which mode dynamics is naturally at- 
tenuated by microscopic effects. One can however imag- 
ine that modes with wave vectors |k| > Dmax/'^ are 
initially small and heavily damped; if this is true they 
will be minimally excited during the dynamics. Under 
these circumstance the condition on tb for an inacces- 
sible D > Dmax ceases to exist (since — )■ oo in this 
case) and we are left only with the condition on the group 
velocity and linear band width, which is similar to what 
is discussed in previous investigations. For D < Dmax 
Td is finite and physically relevant. 
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AFOSR, USA, under the grant FA9550-06-1-0345. 
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